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Spherically symmetric problem on the brane 
and galactic rotation curves 

Alexander Viznyulcl and Yuri ShtanoMl 
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We investigate the braneworld model with induced gravity to clarify the role of the 
■ cross-over length scale I in the possible explanation of the dark- matter phenomenon 

o: 

in astrophysics and in cosmology. Observations of the 21 cm line from neutral hy- 

(N 

drogen clouds in spiral galaxies reveal that the rotational velocities remain nearly 

D ■ _„ , 

constant at a value v c ~ 10 —10 in the units of the speed of light in the region of 



the galactic halo. Using the smallness of v c , we develop a perturbative scheme for re- 



constructing the metric in a galactic halo. In the leading order of expansion in v c , at 
the distances r > v c £, our result reproduces that obtained in the Randall-Sundrum 



braneworld model. This inequality is satisfied in a real spiral galaxy such as our 
Milky Way for distances r ~ 3 kpc, at which the rotational velocity curve becomes 



> 

■^j- | flat, v c ~ 7 x 10~ 4 , if £ ^$ 2 Mpc. The gravitational situation in this case can be ap- 

proximately described by the Einstein equations with the so-called Weyl fluid playing 



the role of dark matter. In the region near the gravitating body, we derive a closed 



o 

system of equations for static spherically symmetric situation under the approxima- 



tion of zero anisotropic stress of the Weyl fluid. We find the Schwarzschild metric 



to be an approximate vacuum solution of these equations at distances r < \Jr g l 2 . 
The value I ^$ 2 Mpc complies well with the solar-system tests. At the same time, 
in cosmology, a low-density braneworld with t of this order of magnitude can mimic 
the expansion properties of the high-density LCDM (A + cold dark matter) universe 
at late times. Combined observations of galactic rotation curves and gravitational 
lensing can possibly discriminate between the higher-dimensional effects and dark 
matter. 
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I. INTRODUCTION 

The idea that our familiar four-dimensional space-time is a hypersurface (brane) of a 
five-dimensional space-time (bulk) was under detailed elaboration during the last 

decade. According to this braneworld scenario, all matter and gauge interactions reside 
on the brane, while gravity can propagate in the whole five-dimensional space-time. An 
observer in this scenario is in direct contact only with the induced metric on the brane. 

A ponderable argument in favor of the braneworld theory is the ability to solve outstand- 
ing problems of modern cosmology and astrophysics. The generally accepted LCDM (A + 
cold dark matter) model is in a good agreement with most of experimental data; however, 
up to now, no non-gravitational evidence for dark matter has been reliably found. It was 

nfinQ n 

argued [4j, la, |6|, [7j (see also |8| for a somewhat different approach) that a modified theory of 
gravity based on the Randall-Sundrum braneworld scenario (jj can explain the observations 
of the galactic rotation curves, while observations of the gravitational lensing can possibly 
discriminate between the higher-dimensional effects and dark matter. In a recent work 9(, 
the authors argued that the virial-theorem mass discrepancy in clusters of galaxies can also 
be accounted for in the Randall-Sundrum model. However, this version of braneworld model 



cannot explain the dark-matter phenomenon on the cosmologica 



scale without introducing 

other components in the theory. At the same time, as shown in [10j |. a low-density braneworld 
with induced gravity can mimic the expansion properties of the LCDM model. In particular, 
a universe consisting solely of baryons with f2b — 0.04 can mimic the LCDM cosmology with 
a much larger 'effective' value of the matter density Q m ~ 0.2—0.3. This effect becomes pos- 
sible owing to the presence of a fundamental cross-over length scale I in braneworld models 
with induced gravity, which is absent in the Randall-Sundrum braneworld model. 

The aim of this paper is to investigate the properties of the braneworld model with induced 
gravity to clarify further the role of the scale £ in the possible explanation of the dark- 
matter phenomenon in astrophysics as well as in cosmology. From the general consideration 



10, LL1[ , one expects the braneworld model to resemble the general relativity at the distances 
t ~ y r g £ 2 , thus giving an approximate Schwarzschild metric in vacuum with r g in a role of 
the Schwarzschild radius. At large distances, the modified gravitational field equations can 
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be used to explain the properties of the galactic rotation curves. 

Observations of the 21 cm line from neutral hydrogen clouds in spiral galaxies reveal that 
the rotational velocities remain nearly constant in the halo region at values v c ~ 10~ 3 — 
in the units of the speed of light. Starting from this observational fact, the exact metric in 
a galactic halo was reconstructed in frames of the Randall-Sundrum braneworld model in 
|6|, |7|, which was used to account for this phenomenon by the effects of higher-dimensional 
gravity without dark matter. We extend the results obtained in these works to the case of 
braneworld model with induced gravity. Using the smallness of v c , we develop a perturbative 
scheme for reconstructing metric in a galactic halo and demonstrate that, in the leading 
order of expansion, at the distances r > v c £, our result reproduces that obtained for the 



Randall-Sundrum braneworld model 



in y,y. 



The gravitational situation in this case can 
be approximately described by the Einstein equations with the traceless projection of the 
bulk Weyl tensor to the brane playing the role of an effective stress-energy tensor of galactic 
fluid. To give a prescription for verifying the scenario based on the braneworld models, we 



apply the results of [121 ] . in which it was proposed to use the existing data on rotation curves 
and lensing measurements to constrain the equation of state of the galactic fluid. 

To solve the braneworld equations at small distances to the gravitating source, one should 
impose some additional conditions. As first noted in [13], the nonlocality and nonclosure of 
the braneworld equations are connected with the projection C a b of the bulk Weyl tensor to 
the brane. It, therefore, seems reasonable to impose certain restrictions on this tensor in 
order to obtain a closed system of equations on the brane. In this way, several classes of 
vacuum static spherically symmetric solutions on the brane were obtained for the Randall- 
Sundrum braneworld model Q, [l^. Solutions obtained by imposing direct restrictions on 
the form of the metric (such as g tt = ~9rr Ql) a ^ so correspond to some evident conditions 
on C a b. Finding the correct boundary condition of this kind, one can turn the braneworld 
model into a complete nonlinear local theory of gravity viable in all physical circumstances. 1 

The simplest way of restricting the tensor C a b consists in setting to zero its (appropri- 
ately defined) anisotropic stress. 2 This condition is fully compatible with all equations of 



1 Another possibilities to restrict the space of solutions of braneworld equations in a spherically symmetric 
case was discussed in 



15 



17 



2 Earlier attempts to simplify the braneworld equations by entirely neglecting the contribution from C a b 
proved to be incorrect since this condition is incompatible with the equation that follows from the Bianchi 



4 



the theory and leads to a brane universe described by a modified theory of gravity with an 
additional invisible component (Weyl fluid) having nontrivial dynamics. A one-parameter 
family of boundary conditions (with zero anisotropic stress as a particular case) was success- 
fully employed in l_8j to derive an exact system of equations describing scalar cosmological 
perturbations on a generic braneworld with induced gravity. 

In the present paper, we derive a closed system of equations for static 3 spherically sym- 
metric situation in the case of zero anisotropic stress in the tensor C ab (the minimal boundary 
condition in the terminology of 18|). We find the Schwarzschild metric to be an approx- 



imate vacuum solution of these equations at distances r < ^/r g £ 2 , which justifies using of 
the minimal boundary condition in the vicinity of gravitating objects. 

The paper is organized as follows. In the next section, we remind the field equations of 
the braneworld model with induced gravity. In Sec. IHH we discuss the generic properties of 
the braneworld equations. In Sec. IIV} we describe the cosmic mimicry model. In Sec. |Vj we 
reconstruct the metric in a galactic halo with flat rotation curves. In Sec. IVIl we derive and 
discuss the system of equations for the static spherically symmetric case with the minimal 
conditions for the Weyl fluid. In Sec. IVHl we find vacuum solutions in the weak field limit. 
In the last section, we present our conclusion. 



II. EFFECTIVE FIELD EQUATIONS 



We start from the simplest generic braneworld model described by the action 



S = M i 



(ft-2A)-2 / /V 

bulk J brane 



/ K + I (m 2 R -2a)+ I L 

J brane J J brane J brane 



(h ab , <p) 



Here, 1Z is the scalar curvature of the metric g a b in the five-dimensional bulk, and R is the 
scalar curvature of the induced metric = — n a n b on the brane (n a is the vector field 
of the inner unit normal to the brane) which is assumed to be a boundary of the bulk space, 
and the notation and conventions of 20[] are used. The quantity K = h ab K a }y is the trace of 



identity [see Eq. ([5]) below] . It was first noted in 13| in the context of the Randall-Sundrum model that 
neglecting C a b in the spherically symmetric situation leads to a severe restriction on the energy-momentum 
tensor of perfect fluid, namely, it only admits solutions with constant matter density (incompressible fluid). 
Below, we demonstrate that this result is valid also in our more general braneworld model with induced 
gravity. 

For a discussion of the gravitational collapse on the brane see [19( . 
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the symmetric tensor of extrinsic curvature K ab of the brane. The symbol L(h ab , 0) denotes 
the Lagrangian density of the four-dimensional matter fields <ft whose dynamics is restricted 
to the brane so that they interact only with the induced metric h ab - All integrations over 
the bulk and brane are taken with the corresponding natural volume elements. The symbols 
M and m denote the five- dimensional and four- dimensional Planck masses, respectively, A 
is the bulk cosmo logical constant, and a is the brane tension. 

Action ([I]) leads to the Einstein equation with cosmological constant in the bulk: 

Gab + Agab = , (2) 

while the field equation on the brane is 

m 2 G ab + ah ab = T ab + M 3 (K ab - h ab K) , (3) 

where T ab is the stress-energy tensor on the brane stemming from the last term in ([1]). The 
Einstein equation in the bulk ([2]) and the Gauss-Codazzi identities imply the conservation 
of the stress-energy tensor of matter on the brane, namely, V a T ab = (see, for example, 
2ll|). Here, V a denotes the covariant derivative on the brane associated with the induced 
metric h ab . 

By using the Gauss-Codazzi identities and projecting the field equations to the brane, 
one can obtain the following effective equation lOj: 



where 



are convenient parameters of the braneworld model, 

A a 2 

ArS= 2 + 3M^ (6) 
is the value of the effective cosmological constant in the Randall-Sundrum model, 

Qab = ~EE ab — E ac E c b + - ^E cd E cd — ^E 2 ^j h ab (7) 

is a quadratic expression with respect to the 'bare' Einstein equation E ab = m 2 G ab — T ab on 
the brane, and E = h ab E ab . The symmetric traceless tensor C ab = n c n d C acb d is a projection 
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of the bulk Weyl tensor C a bcd which carries information about the gravitational field outside 
the brane. The tensor C ab is not freely specifiable on the brane, but is related to the tensor 
Q a b through the conservation equation 

V a {Q ab -M e C ab ) = 0, (8) 

which is a consequence of the Bianchi identity applied to (j3J). 

The following famous braneworld theories are related to important subclasses of action 
©: 

1. The Randall-Sundrum model is obtained after setting m = in (JT]). The corre- 



sponding effective equation on the brane was derived in 131 ] : 



where 



G a b + ^Rshab — 3 ^ ab ^ ab ~~ ^ ab ' (^) 



Sab = \TT ab - T ac T% + l - (r cd T cd - l -T 2 ) h ab . (10) 



In this model, only the right-hand side of the Einstein equation is modified. 

2. The Dvali-Gabadadze-Porrati (DGP) model [s| corresponds to the special case where 
both the cosmological constant in the bulk and the brane tension vanish, i.e., A = 
and a = in (TjQ). The corresponding effective equation 

Gab = J^Qab — Cab (H) 

does not contain a linear contribution from the stress-energy tensor T a b, which makes 
a perturbative analysis in this theory problematic. 

3. Finally, general relativity, leading to the LCDM cosmological model, is formally ob- 
tained after setting M = in (pQ). In this case, the effective equation Q a b = has an 
obvious solution 

G a b = — 9 T a t) . (12) 
III. GENERIC PROPERTIES OF BRANEWORLD GRAVITY 



The theory described by action ([T]) has two important scales, defined in namely, the 
cross-over length scale £, which describes the interplay between the bulk and brane gravity, 
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and the mass (energy) scale k, which determines the role of the brane tension in the dynamics 
of the brane. Thus, we can expect the properties of braneworld gravity to be different at 
different scales. 

Indeed, looking at (J3|, we notice that expression for Q a b is quadratic in the curvature 
as well as in the stress-energy tensor. Therefore, in the region of high matter density and 
curvature, the tensor Q ab dominates in (j2J), and the gravitational law is approximated by 
the 'bare' Einstein equations m 2 G a b — T ab = 0. And, vice versa, the contribution from Q ab 
in (jlj) is insignificant on sufficiently large length scales, where the curvature is small, and 
the braneworld theory on those scales should again be approximated by the Einstein gravity 
with different gravitational constant. The tensor C a b in this case plays the role of some 
effective stress-energy tensor. 

To give some qualitative numerical estimates, we repeat here the reasoning of [h]]. Con- 
sider the trace of the effective equation (J4]): 



where the left-hand side contains terms which are linear in the curvature and in the stress- 
energy tensor while the right-hand side contains the quadratic term Q = h ab Q ab . This 
equation is closed on the brane in the sense that it does not contain the contribution from 
the traceless tensor C a b- 

Suppose that we are interested in the behaviour of gravity in the neighborhood of a 
spherically symmetric object with density p s , total mass Ai s , and radius r s . For simplicity, 
we assume that one can neglect the tensor C a b and the effective cosmological constant in the 
neighborhood of the source. As regards the effective cosmological constant, this assumption 
is natural. Concerning the neglect of the tensor C a b, we demonstrate the validity of this 
approximation in Sec. I VIII by considering the vacuum spherically symmetric situation in the 
braneworld model with C ab restricted by a simple condition of vanishing anisotropic stress. 

Within the source itself, we have two qualitatively different options: an approximate 
solution can be sought either neglecting the quadratic part or linear part of (jlj) and ({TBI . 
We should choose the option that gives smaller error of approximation in ( fl3|) . In the first 
case, neglecting the quadratic part and the effective cosmological constant, we have 




(13) 




(14) 
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In the second case, we neglect the linear part, so that 

Q ab ^0 =► E ab ^0 =► R+ (^-) — 2 T~ P° (15) 

The final expression on the right-hand side of (fl5l) is smaller than the corresponding expres- 
sion in (THj) if 

p a > + 1) 2 — => rg<rg~ ' 2 , 16 
m z [p + l) z m 2 

where we used the relation M. s ~ p s r;?. Thus, we can expect that, in the neighborhood of 

the source, on distances smaller than r* given by ([TBI) , the solution is determined mainly 

by the quadratic part Q a b in ©, which means that it respects the 'bare' Einstein equation 

m 2 Gab — T a b = to a high precision. This effect is sometimes described as the 'gravity filter' 

of the DGP model [3], which screens the scalar graviton in the neighborhood of the source 

making the gravity effectively Einsteinian. Some aspects of this interesting phenomenon 

were discussed in [22]. 

Expression (f!6l) generalizes the length scale [Uj of the DGP model, below which nonlinear 
effects become important, to the case of nonzero brane tension (nonzero (3) and bulk cos- 
mological constant satisfying the Randall-Sundrum constraint Ars — 0. In order to comply 
with observations in the solar system, the value of i should be chosen sufficiently large. In 
the next section, we will see that modeling the homogeneous dark matter phenomenon by 
cosmic mimicry requires (3 ~ —5/4, and in Sec. |V]we will show that matching the galactic 
rotation curves can be realized if i ^ 2 Mpc. For these parameter values, the value of r* for 
the Sun is of the order of 1 pc, which is quite large so that the solar-system experiments are 
in accord with the general-relativistic expectations. 



IV. COSMIC MIMICRY ON THE BRANE 

In our paper jlQl ]. we described a braneworld model in which cosmological evolution 
proceeds similarly to that of the Friedmannian cosmology but with different values of the 
effective matter parameter Q m , or, equivalently, with different values of the effective gravi- 
tational constant G^, at different cosmological epochs. 

Specifically, for a spatially flat universe with zero background dark radiation (C = 0), 
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the cosmological equation stemming from ([I]) can be written as follows: 

2 

/ i \ 2 

A „ _ / / A 1 \ 1 

r2 



6 



\ 3m 2 



+ 



A 
6" 




(17) 



where po an d Hq are the energy density of matter and Hubble parameter, respectively, 
at the present moment of time. The two signs in (1171) correspond to two complementary 
possibilities for embedding the brane in the higher- dimensional (Schwarzschild-AdS) bulk. 



In the effective mimicry scenario 



lOj , the parameters Hq — A/6 and 1 ji 2 are assumed to 



be of the same order, and much larger than the present matter-density term p /3m 2 . The 
mimicry model has two regimes: one in the deep past (where the matter density was high) 
and another during the late-time evolution (where the matter density is low). In the deep 
past, we have 

1 



P~ Po 
3m 2 



> 



6 



and the universe expands in a Friedmannian way 

P 



H 

During the late-time evolution, we have 

P~ Po 



3m 2 



(18) 



(19) 



3m 2 



< 



and the expansion law is approximated by 

H 2 » H 2 + 

where 

a = 



a P~ Po 
a =F 1 3m 2 



m- 



A 
"6 



(20) 



(21) 



(22) 



is the parameter introduced in [10[]. In the case of the cosmological branch with upper sign, 
it is assumed that the coefficient a/ (a — 1) in ( 12TI) and in similar expressions is always 
positive, i.e., a is assumed to be greater than unity in this case. 

One can interpret the result (J2"TT) either as a renormalization of the effective gravitational 
constant relative to its value in the deep past or as a renormalization of the effective density 
parameter: 

,L( 'L>M L('L)i\J 

p LCDM = ^_ p _ (23) 



H 2 « 



LCDM 

Po 



3m 2 



a =F 1 
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Remarkably, the behaviour of the Hubble parameter on the brane practically coincides 
with that in LCDM at low densities. This property was called 'cosmic mimicry' in |l0| . 
A consequence of this is the fact that a low/high density braneworld consisting entirely of 
baryonic matter with ~ 0.04 could easily masquerade as LCDM with a moderate value 
qlcdm 0.2—0.3. Thus, with the upper sign in ( 1231) . we can explain the phenomenon of 
homogeneous dark matter if suppose a ~ 5/4. In this case, taking into account that the 
value of I lies well below the Hubble scale, from (1221) we obtain 

|A|£ 2 ^6a 2 ^9. (24) 

The value of A should be negative for the expression under the square root in (122]) to be 



positive. As was shown in [10|, (3 is also negative in this case, and \(3\ ~ a. The upper 
estimate i < 2 Mpc is obtained in the next section by considering galactic rotation curves. 
The range of redshifts over which this cosmic mimicry occurs is given by < z <C z m , where 
z m in this case can be estimated to be z m ^ 170. 

V. RECONSTRUCTION OF THE METRIC IN A GALACTIC HALO 

In the previous section, we described how the phenomenon of homogeneous dark matter 
can be explained by cosmic mimicry in the braneworlds models with induced gravity. The 
question then arises whether the phenomenon of dark matter can also be accounted for in 
the inhomogeneous situations, in particular, on the scales of galaxies. This is the subject of 
the present section. 

The problem of reconstructing the static metric and dark-matter halos in general relativity 
from the behavior of galactic rotation curves was investigated recently in [l^, 23I 24]. A 
similar procedure was applied to the Randall-Sundrum braneworld model in 5_|, |g, |7| , with a 
conclusion that it can explain the observations of the galactic rotation curves without dark 
matter. The Randall-Sundrum model, however, does not allow for the mimicry property 
and thus cannot explain the homogeneous dark-matter phenomenon. Hence, it is necessary 
to address this issue in frames of the induced-gravity model. 

The traceless tensor C a b drops out completely from the trace of (T4|), which is the only 
closed equation on the brane in the absence of any additional constraints on the tensor 
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251 ]. In the vacuum, the trace of (@J reads 

1 



i • >')/.' • y 2 { R * bRab - s R2 1 1 A 7 = 



(25) 



One can use this equation to reconstruct the metric in the neighborhood of a galaxy by 
assuming approximate spherical symmetry and taking into account the qualitative behavior 
of rotation curves. In the Randall-Sundrum model [which corresponds to setting £ = and 

(3 = in (1251) ]. this procedure was done in aaa. 

The spherically symmetric metric, by which we approximate the situation in a galactic 
halo, has the form 



ds 2 = -f{r)dt 2 + g-\r)dr 2 + r 2 (dB 2 + sin 2 6dip 2 ) . 



In the region of the flat rotation curve, we have the condition 



<j>'(r) 



v- 



(26) 



(27) 



where v c is the rotation velocity, and <f>(r) = ~ log/(r). It can be integrated to give the g o 
component of the metric: 

/W=(-) , (28) 

where r c is the constant of integration. Substituting metric fl2"B"|) with the function <p(r) given 
by fl28l) into equation fT25|) . we obtain the following first-order nonlinear differential equation 
for the function g(r) : 



(rg'-2gf + Arg' 



12(1 + I3)r 2 \ ( 6(1 + /3)r r < 

-%g 1 + — To — 



+ 4 



1 + 



12r 2 (l + ( 5-2A RS r 2 



-2r 2 Q/) 2 - 24rgg> + 8g 2 + 8rg' (l - 3(1 - 8g (l 



6(1 + py 
I 2 



{g') 2 -12rgg' + 16g 1 



3(1 + py 



+ 4v b c g(rg'-4g)+4v 8 c g 2 = 0. (29) 



This equation looks rather complicated. However, we can take into account that the 
quantity v 2 ~ 10 -7 is very small for typical galaxies and thus represents a convenient 
parameter for asymptotic expansion. Therefore, we can look for the solution of (1291) in the 
form of expansion in powers of v 2 : 



g{r) = go{r) + v 2 c gx{r) + v%(r) + ... 



(30) 



12 



where the functions g , g\, g 2 , ■ ■ ., should be determined step by step by solving (129]) with 
accuracy corresponding to the order of v c at every step. First, we should find go( r ) ; which 
is the solution of ( 1291) with all powers of v c neglected. Introducing %{ r ) = r 9o( r ), we can 
rewrite the equation for go(r) in the following form: 

f (x' ~ y + 2) 2 + 48r 2 [(1 + /?)(1 - X ') - 2A RS r 2 ] = . (31) 

In the Randall-Sundrum limit, which implies I = and /? = 0, we should have Xi( r ) = 
r + a — 2ARsr 3 /3 as a solution of (13*11) . where a is some arbitrary constant of integration. 
On the other hand, the function Xzi 7 ") — r + ^ r3 > with another arbitrary constant b, sets 
expression in the first parentheses to zero. One can easily observe that the function x( r ) = 
r — 2Ars^ 3 /3(1 + (3) is an exact particular solution of (13T|) . setting to zero both components 
on the left-hand side separately. Thus, for the zero-order approximation, we have 

*M = l-^. (32) 

The expression 2Ars/(1 + f3) here represents an analog of the cosmological constant in 
general relativity. At the scales at which galactic rotation curves are observed, the contri- 
bution from the cosmological constant should be negligibly small, namely 

thus allowing us to neglect it in the future. We should keep it for a moment, to estimate 
corrections in the second order of v c . 

Using (J3"0|) and ("3~2"]) . we find the differential equation for gi(r) : 



^ A RS £ 2 [ 1 4A RS r 2 ^ 3(1 + /?) 2 

rgx + gi = -1 + 



3(1 + /?)* 

which can easily be integrated to give 
^i(r) = - 1 + — 



1 - 



02 



r J 3(1 +/3)' 



3(1 + /?) V A RS £2 



4A RS r 2 / 3(1 + PY 



9(1 + /?) V A 



(34) 



(35) 



where C\ is some constant of integration having the dimension of length. 

In the cosmic-mimicry theory which, as described in the previous section, explains the 
dark-matter phenomenon on the cosmological scale, the value of £ lies well below the cos- 
mological length scale. Thus, we have 
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Conditions (|33|) and (|36|) allow us to write the function g(r) at this stage in the form 

2A RS r 2 2 / CA 

^) wl -3crT^" Ve l 1 + Tj ■ (37) 

It is important to note that, under conditions (133]) and (|36|) . the difference between 
g(r) in our model and the same quantity in the Randall-Sundrum braneworld model in this 
approximation in v 2 c is only in the appearance of the factor (1+/?), which simply renormalizes 
the cosmological constant. 

To find corrections of the forth order in v c , we shall completely neglect the contribution 
from the term containing A^g. This can be done under the condition 

(38) 

which, obviously, is more restrictive than (133|) . but still is satisfied for a galaxy The result 



is 



g(r) & 1 — v 



2 v „ , r 



,2 



2V1 



Pvt / d C 2 



4(1 + /?)r 2 V 2r 4r 2 



+ , (39) 



in which the constant of integration Ci is renormalized as compared with the previous 
approximation. 

Although the last term in (l3T?j) is proportional to t; 4 , at too small distances it can turn 
out to be larger than v 2 c . This will obviously destroy our perturbative analysis. To avoid 
such situation, we restrict ourselves to the region 4 r 2 3> £ 2 v 2 /|l + (3\. 

Leaving only corrections of the leading order in (|39|) . we can write 

g{r)Ml-vl(l + Q\ , (40) 



r 

which is a general reconstruction of the metric in our braneworld model in the region 

\rw\ <<r<< ^r (41) 

with condition (1281) . 

In order to satisfy the left inequality in (1411) in a real spiral galaxy such as our Milky Way 
with the distance r ~ 3 kpc at which the rotation velocity becomes constant of magnitude 



For simplicity, we assume the length scale \C±\ to satisfy the condition \C\\ < v c t/ + (3\, so that 
|Ci| /r is small. However, we would like to note here an interesting possibility that a large constant \Ci\ 
determines another interesting scale in a galactic halo. 
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v c ~ 7 x 10 -4 , we should have £ < y|T+~/3j -U" 1 r ~ 2 Mpc, where we took into account 
the mimicry estimate (3 = —a ~ — 5/4 (see the end of Sec. UVp . In this paper, we restrict 
our investigation to this case. For magnitudes of the cross-over length scale I essentially 
larger than this value, one cannot use the perturbation expansion in powers of ir?, and it is 
necessary to solve the nonlinear equation (1291) . 

We can substitute ff28|) . fT40l) into (jlj to compute the tensor C a b- The result is 

C,»fi±M, C7 r ,(i±M(%-l), ., = > yW _ (l±» | (42) 



where we have neglected the contribution from terms of order u 4 , which is correct in region 
(jUJ). Again, we observe that, in the leading order of expansion, our result reproduces that 
obtained for the Randall-Sundrum braneworld model in 6|,|7j]. This is what we have expected 
from general consideration of Sec. IIII1 on sufficiently large length scales, the contribution 
from the quadratic expression Q a b in (j4j) is negligibly small. The gravitational situation 
in this case can be approximately described by the Einstein equations with the tensor C a b 
playing the role of some effective stress-energy tensor T ab : 

8nG^T ab = — - — p C ab = G ab , (43) 

where Gn is the gravitational constant adopted by an observer on these scales. 

Using fT42l and (j4~3]) . we determine the effective energy density p(r) and radial and trans- 
verse pressures p r (r) and pt(r), respectively, inside the galactic halo: 

8vrG N p(r) « V l , 8vrG N p r (r) « V l (l - ^ , 8nG N p t (r) « ^ . (44) 

We note that p r (r) 7^ pt{r), which indicates the presence of a significant anisotropic stress 
in the tensor C a b [see definition ( 14T1) below] on galactic scales. 

It was argued in l_2j that combined observations of galactic rotation curves and gravi- 
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tational lensing can yield the profile of pressure inside a galactic halo. The authors of 
proposed to use the pseudo-masses m RC (r) and m lens (r), which can be obtained for the same 
galactic halo from the treatment of rotation curves and gravitational lensing in the context 
of dark-matter paradigm (which implies negligible pressure), to define the x-I&qXoy which 
quantifies the deviation from the predictions of the CDM model: 

Y[ W (r)l = mLs{r) = 2 + Mr) (45) 
m^ c (r) 2 + 6a;(r)' 1 ' 
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where 



p r (r) + 2p t (r) 
3p(r) 



(46) 



For the braneworld model, using (jUj), one obtains w«l/3 and x « 3/4, meaning that 
deflection angles, which can be computed in our model, will be around 75% of the usual 
predictions based on dark matter. This result coincides with the analysis made in s|. 



VI. MINIMAL CONDITION FOR THE WEYL FLUID ON SMALL SCALES 

In the previous section, we have shown that the contribution of Q a b to (jl]) is negligibly 
small on sufficiently large spatial scales, which makes it possible to explain the flatness 
of galactic rotation curves as an effect of higher- dimensional gravity. On such scales, the 
braneworld theory reduces to the Einstein gravity with the tensor C a b playing the role of 
some effective stress-energy tensor. 

To demonstrate that the contribution from C a b to the gravitational dynamics is insignif- 
icant at small spatial scales, we should solve Eq. (j4j) at small distances from the source of 
gravity. The problem arising here is that equation (j3J) is not closed on the brane in the sense 
that the dynamics of the symmetric traceless tensor C a b on the brane is not determined 
by the dynamics of matter alone. An approximation usually taken by several authors to 
overcome this problem consists in imposing certain conditions on the tensor C a b directly on 
the brane so as to close Eq. (T4]). Such conditions on the tensor C a b should be compatible 
with the conservation equation (J5J) and should also leave the braneworld theory compatible 
with Einstein's general relativity in a wide range of physical situations. 

In general, the tensor C a b can be decomposed through an arbitrary normalized timelike 



vector field u a on the brane (see 



2d) = 



m 2 C ab = y {h ab + Au a u b ) + 2v c (a u b) + it c ab . (47) 

Here, the covector v c a and traceless symmetric tensor i^ c ah are both orthogonal to u a . The 
tensor C a b in this case is regarded as the stress-energy tensor of an ideal fluid with equation 
of state like that for radiation but with nontrivial dynamics described by Eq. OH]). In the 



literature, it is called 'Weyl fluid' 27| and, in the cosmological context, 'dark radiation' 28]. 
The quantity pc has the meaning of its density, v c a is its momentum transfer, and 7r^ b is its 
anisotropic stress, as measured by observers following the world lines tangent to u a . The 
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stress-energy of this ideal fluid is not conserved due to the presence of the source term Q ab 
in ©. 

Equation (jHJ) gives evolution equations for the components pc and v$. However, there 
are no evolution equations for the tensor fields 7Q, on the brane, which is a manifestation 
of the nonlocality of the physical situation from the braneworld viewpoint. Thus, boundary 
conditions for the brane-bulk system can be specified by imposing additional conditions on 
the tensor C ab . 

As the first condition on the tensor field C ab , we demand that it has a normalized timelike 
eigenvector field Uq, so that C a b u b c oc u%. Then this vector field can be used as vector u a in 
decomposition (jH|), which was not specified up to now: 

m 2 C ab = ^ (h ab + Au c a u c b ) + v c ab . (48) 

Note that the covector component v c a vanishes in this case because of the eigenvector property 
of Uq. The Weyl fluid is now described by the timelike eigenvector field Uq, which can be 
interpreted as its four- velocity, and by the quantity pc, which has a meaning of its rest-frame 
density. 

Now, in solving the spherically symmetric problem in the neighborhood of the source, we 
assume that the anisotropic stress can be neglected in the natural decomposition (148!) : 

<t = • (49) 



This is a minimal boundary condition for the brane-bulk system in the terminology of 18]. 
Under condition fj49|) . the normalized eigenvector field Uq is unique modulo orientation. 

Equation ([8]) gives the evolution equations for the components pc and Uq. Hence, fl3j), 
together with equations for material fields, leads to a complete set of equations. 

In what follows, we assume (T49l) to be the case at small distances to the source of gravity 
and solve the vacuum static spherically symmetric problem with the additional condition 
( pf9l . On large scales, this conditions is violated [see Eq. (1441) ]. 

A general static spherically symmetric metric on the brane has the form 

ds 2 = -f{r)dt 2 + h{r)dr 2 + r 2 {dd 2 + sin 2 Bdy 2 ) . (50) 
It leads to the following nonzero components of the Einstein tensor: 
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G r 



rfh 



G% = 



f" f'h' 



r + r 



h! 



(52) 



(53) 



2fh 4fh 2 4f 2 h 2rfh 2rh 2 ' 
The interior of static objects (such as stars) is naturally described by the stress-energy 
tensor of perfect fluid T a b = pu a Ub + p(h a b + u a Ub). To be compatible with the space-time 
symmetry, the fluid four-velocity u a must be aligned with the static Killing vector field £ a : 



u 



(54) 



so that matter is described by the two functions p(r) and p(r), and the coordinate compo- 
nents of its stress-energy tensor are given by 

~l>\ I i . u ^ 

Ta 



( 



(55) 



p{r), 

, p(r)^ J 

The spherical symmetry and static property imply that the Weyl fluid four-velocity is 
also aligned with the static Killing vector field: 5 

d 



(56) 



so that the Weyl fluid is described by the single function pc{r): 

\ 



m 2 C a 
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-Pc{r) 







(57) 



One should note that C r r = C e e as a result of our boundary condition f j49l) which requires 
vanishing of the anisotropic stress. The tensor C a b manifests itself directly as dark radiation 
in the Randall-Sundrum model [see Eq. (jHJ)]; however, in a braneworld model with induced 
gravity, which is described by the effective equation (j3J), its influence on the gravitational 
dynamics is more complicated, and new effects can be expected. 
To proceed further, we rewrite Eq. (j3J) as 



n^ofc + C a b + A RS h a b 



(58) 



5 The four- velocity of the Weyl fluid certainly will play an important role in a non-static situation. 
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Here, Q a b is the quadratic expression with respect to the nonzero diagonal components of 
E ab given by ([7]). Introducing 



x = -E\ , y = E r 



z = E y 



and 



7 



3M 6 (1 + (3) 



7l = 3M 6 



P + Pc 



Ars , 72 = 3M 6 



P + \pc 



A 



we obtain from (1581): 



x 2 — y 2 — z 2 + 2yz = 7X + 71 , 
x 2 - y 2 + z 2 + 2x2 = 7y + 72 , 
x 2 + y 2 + xy + xz — yz = jz + 7 2 , 



(59) 

(60) 

(61) 
(62) 
(63) 



which is a system of algebraic equations with the components of E a p playing the role of 
unknown variables. 

Because of the condition C r r = C e g, the same constant 7 2 appears on the right-hand 
side of (1621) and (1631) . This allows us to solve this system of algebraic equations in the most 
general case. Subtracting (1621) from (|63|) . we obtain 



(y-z)(x + 2y + z + 1 ) = 0, 

which admits two alternative solutions: y = z or x + 2y + z + / y = 0. 
First, consider the case y = z. The solution of (IBTj) . f[6"2"l) is then: 



(64) 



6(1 + /?) 

£ 2 



p + pc 3(1 +/3) 2 

H ^ A RS 



6(1 1 



£ 2 



p-p+2p c /3 , 3(l+/3) 2 _ a 

2m 2 i m ^ 



i-RS 



P+PC , 3(l+/3) 2 A 
"772 1 72 ^RS 



r e | 6(l + /3) _ 1 6 
w H ^ — — + -7 



p-p+2p c /3 , 3(l+/3) 2 _ a 
2m 2 ^ P 1 



RS 



77T 



P+PC I 3(l+/3) 2 _ a 

™2 ~r »2 



(65) 



(66) 



(67) 



^RS 

The second solution of (164)) . namely, x + 2y + 2 + 7 = 0, leads to the result 
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6(1 + (3) _ J_ p ± 2 /p + pc/3 3(1 + (3f 
m 2 i V m? 



GT r + = ~ 2 P ± 7 V ^ + Ars , (69) 



p-p+2p c /3 , 3(l+/3) 2 a 

— w 1 p A RS 



e 6(1 + /?) 1 2 
e H « = — ~p ± - 

£ / p+Pc/3 _ 3(l+/3) 2 , A 

,2 (2 T ^i-RS 



(70) 

It is remarkable that the effective braneworld equations in a static spherically symmetric 
situation in our model have the form of Einstein equations with a special modification on 
the right-hand side. We should also note that the system of equations ( |65l) - (|67j) matches 
with the two branches of cosmological equations ( 1171) of the braneworld model derived in 



29j. This means that writing the components G a p for the Friedmann-Robertson- Walker 
metric and substituting them to the left-hand side of (!65j) - (l67|) . one obtains the cosmological 
equations ffTTl) . As for the branches described by equations (168]) - (ITU]) , their connection with 
cosmology is not clear. Because of this, we restrict our investigation only to the analysis 
of system (l65 i) -(l67 l) . Following the cosmological nomenclature, we call the branch with the 
lower (' — ') sign Branel, and the branch with the upper ('+') sign Brane2. These two 
branches have a geometric origin: t hey correspond to the two possible branches of solutions 
in the five-dimensional bulk space [291 ] . 

It is convenient to introduce new functions 



and a new constant 



g = p±2fo( A /l + £±^-l) , (71) 

■y _|_ p-P+2pc/3 

V = V TltPL | . 2PL -1| , (72) 

<l _|_ P+Pc 



A=ifl+/3±i) , (73) 



where pl and £ are constant parameters defined by 

3m 2 L 



PL = ^T> £ = y ( 74 ) 



and L is a new length scale 



L 2 ~ £ 2 + 3m 2 6 ' ^ } 

We assume the expression on the right-hand side of (175]) to be positive to avoid singularity 
in the square root of (ED) , (|72l . hence, also in system ( |65l) - (|67l) . for positive values of p and 
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Pc- 6 To avoid a similar singularity in (|68|) - (!70|) . one would require the opposite condition, 
namely l/£ 2 + a/3m 2 — A/6 < 0. Thus, one can think that the choice between the two 
systems of equations (l65l) - (|671) and ( l68l) - (170l) is determined by the value of the expression 
on the right-hand side of ( J75i) . If this expression is equal to zero, then these two systems 
coincide. 

The signs of o (or (3) and A are not yet specified, but one should note that the above 
condition leads to the constraint (3 > — 1/2 + A£ 2 /12. Transition to the DGP model in our 
equations is realized by setting a = and A = 0, which implies (3 = and £ = 1. In the 
cosmic-mimicry model with (3 ~ —a ~ —5/4, which we discussed in Sec. ITVl the value of £ 
can be found to be 

1 -. « 4. (76) 
^ |1 + /3| v ; 

In principle, the density of the Weyl fluid pc, unlike that of realistic matter density, is 
not restricted in its sign. However, the above condition of absence of singularity restricts 
the range of boundary conditions for pc( r ) which we shall specify in the next section. 

In the new notation (I7ll-(l73l). system (l65l)-( l67j) is written as follows: 

(77) 

(78) 
(79) 

The quantities q and V play the role of effective energy density and pressure in the usual 
Einstein equations describing the static spherically symmetric situation, and A is the effective 
cosmological constant. The formal solution of (FT?]) - (179]) is well known: 

ds 2 = - e mr) dt 2 + + r 2 , dQ 2 + gin 2 ^ 2N (g0) 

1 — 2p(r)/r v 



-G\ 


= A + 


1 


Kjc r — 


-A + 


9 ' 

m z 


Lx e — 


-A + 


9 



where 



= p(r)+r 3 P(r)/2m 2 -Ar 3 /2 
(?{r} r[r-2fi(r)] ' 1 ' 



6 These singularities directly correspond to the 'quiescent' singularities in braneworld cosmology discussed 
in [io | . 
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The distributions of the densities of usual matter and Weyl fluid are determined by the 
energy-momentum conservation and by the analogue of the Tolman-Oppenheimer-Volkoff 
equation [which, in our case, is a simple consequence of the Bianchi identity applied to (1771)- 
(1791)]. In Einstein's theory, the energy-momentum conservation equation and the Bianchi 
identity give the same result, which is precisely the Tolman-Oppenheimer-Volkoff equation. 
In our model, due to the presence of effective g and V on the right-hand side of (!77j) - (!79j) . 
which do not coincide with the usual p and p, we have two different equations: 

p' + <P'(p + p) = 0, (83) 

I + C) Pc + (1 + C)P' = 8 -y<P' , (84) 

C = P + Pc/3 ' PL . (85) 

P + PC + PL 

For a given equation of state p = p(p) and some boundary conditions for p(r) and pc{r), 
one can integrate (at least, numerically) (18~T|) . (1521) . (I8"3"j) . and (18"4"|) to obtain the metric and 
matter density distribution inside a massive object. As for the boundary conditions, it looks 
natural to specify the values of p and pc at the center of the spherically symmetric gravitating 
body. At the surface of this object, the interior solution matches with the exterior vacuum 
solution of the above equations. 

It should be noted here that Eq. ( 154"1) gives the condition of incompressible fluid p'(r) = 
for the ordinary matter if we set pc{r) = 0, which corresponds to a complete neglect of the 



where 



tensor C a b in the initial field equations 13|, |3jj . Allowing for nonzero density of the Weyl 



fluid replaces this restrictive condition on p(r) by a differential equation determining the 
behavior of pc( r )- 

We see from ( 1771) - (1791 that the effective gravitational constant is given by 87tGn = 1/m 2 . 
For £ w 2 Mpc and £ ~ 4, the value pi in (1741) can be estimated to be pl ~ 10~ 22 g ■ cm -3 . 
This density scale lies far below the density of compact objects, such as stars and planets. 
Thus, inside such objects, we expect only slight modification of Newtonian dynamics due 
to the bulk effects [g{r) ~ p(r) and V{r) w p{r) if pcij-) < 1]. The vacuum situation 
corresponding to this case will be studied in the next section. 

The local disk density of a typical spiral galaxy is pjisk ~ 3—12 x 10~ 24 g • cm -3 , and, 
therefore, the bulk effects inside the disk are significant. The solution of the interior problem 
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for such objects is of special interest, but is beyond the target of this article. The solution 
to this problem will hopefully give a relation between the asymptotic rotational velocity v c 
and the gravitational mass of a spiral galaxy. 



VII. VACUUM SOLUTIONS 



A. General equations and boundary condition 

In the vacuum, which is defined by the conditions p(r) = and p(r) = 0, the effective 
matter density g(r) and pressure V(r) are not necessarily equal to zero due to the possible 
presence of the Weyl fluid component pc( r ) : 

ft^(r) = ±2£p L ( Vl + Pc(r)lp L - l) , (86) 

P„(r> = ^p L ( l±M ^ - l) . (87) 
VV 1 + Pc{r)/pL J 

One should note that the effective equation of state in this case 

V™_ 1 ( 2 _ \ 1 1 T Q v J(2£pl) (g8) 
Qvac 3 V y/T+ p c /p L J 3 1 ± ^ vac /(2^p L ) 

differs from that of radiation. Specifically, V vac / g va , c < 1/3 if pcij) is positive. 

In addition to equations (1771) . (1781 . we can use Eq. instead of f!79|) in the vacuum 
case: 

P -^-^pi c = Ap c <P>. (89) 
Pc + Pl 

The key issue concerns the choice of the boundary condition for pc{r). As we noted above, 
it is natural to specify the value of pc at the center of the spherically symmetric gravitating 
body. However, as the interior problem is rather complicated and because we are not going 
to consider it in this paper, it is reasonable to parameterize the boundary condition by the 
value of pc in the vacuum at some radius R. For a massive body, this could be the radius of 
its surface, from which the vacuum solution starts. Thus pn = pc(-R) is a new free parameter 
of our model, defining the boundary condition for the Eq. (1891) in vacuum. It is restricted 
by the condition pn> ~Pl, which is required for the expression under the square root of 
the right-hand side of fl86|) and flHTI) to be positive. 
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Eq. (!89|) can be solved exactly, relating the energy density of the Weyl fluid pcij) and 
the function <p{r) entering metric fl80l) : 

Pc(r) PR _ e A[4, R -^r)\ ^ (gg) 



[pL + Pc(r)] 2 (pl + Pr) 2 ' 

where 4>r = <f)(R)- If the boundary value p# = 0, then pc( r ) = for r > R. In this case, 
Qvac{ r ) = V va ,c(r) = [see (1861) and (IHTj) ]. and we obtain the usual vacuum Einstein equations 
with a cosmo logical constant A. Thus, our braneworld model admits a Schwarzschild-(A)<iS'4 
space as an exact vacuum solution. The interior counterpart of this exterior solution in the 
context of the Randall-Sundrum braneworld model was analyzed in HQ. 
To solve equations f!77|) and f!78|) in the vacuum, we introduce the function 

P 2<£(r) 

A(r) = , w . (91) 

v ; l-2//(r)/r v ; 

Then one has 

% = X±^(f^-l), (92) 
m \y Pl J 

l -(l- 2 A- = PC (93) 

r V r J A 3m2 ^1 + p c /p L ' 

This nonlinear system in its full generality looks quite complicated. But, in fact, one can note 

that the gravitational field is usually weak in the neighborhood of an astrophysical object. 

So, in analysis of ( |92l) . ( l93l) . we restrict ourselves only by the Newtonian approximation. 



B. Newtonian approximation 

First, we note that the function 0(r) is an analogue of the Newtonian gravitational 
potential in this case. Under the condition 

|0(r)|«l, (94) 
Eq. (|90|) implies an almost constant value of pc{r) : 

Pc(r)^p R . (95) 
In this approximation, Eq. (19"2"|) can be integrated: 
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where r g is the integration constant and is an analogue of the Schwarzschild radius in general 
relativity. 

The Newtonian approximation requires also the condition \fi(r)/r\ <C 1. Assuming it to 
be satisfied, we can solve Eq. ( 1931) : 

A(r) « e ±r2/r N « 1 ± ^ , (97) 

where 

rN "^' \ P rIpl\ m 
determines the radius at which the Newtonian approximation fails. The integration constant 

that arises in solving (|93|) is absorbed by rescaling the time coordinate. The parameter is 

unambiguously defined by and can be taken as a new free parameter of our model. The 

Newtonian potential in this case is defined by 

r Xr 2 2fr 2 [V 1 + Pr/Pl ~ 1 ) 
e 2 ^) » 1 + 20(r) « 1 - ^ - ^ ± ■ V 7 ■ 7 - , (99) 

r 3 L 2 v 71 + Pr/Pl 

and the region of applicability of the above result is 

r g < r < r N . (100) 

One can easily verify that, under condition (11001) . our initial assumptions \4>{r)\ <C 1 and 
|/x(r)/r| 1 are satisfied. If Pr/ pi is of order unity, then r N ~ L/y/£. Equations fl96|) and 
( |99l) then describe the Schwarzschild metric in the region r 3 ^ £,r g £ 2 . For ^ given by (J76l) . 
this condition is compatible with condition ( fl6l) . thus justifying our initial assumption that 
bulk effects is insignificant at sufficiently small scales. 

A remarkable feature of our model in comparison with the Randall-Sundrum braneworld 
model is the appearance of a new length scale I [which goes to zero in the limit m — > 0; 
see ( !75|) ] defining the radius r* ~ r g £ 2 [see (Tl6l) ] up to which the Newtonian gravity works. 
Beyond the radius r*, our model predicts transition to gravity with different properties, 
for which the influence of the projection of the bulk Weyl tensor C a b on the dynamics can 
be significant, thus allowing to explain the flatness of galactic rotation curves, as we have 
demonstrated in Sec. [V] 
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VIII. CONCLUSION 

In this paper, we continued the analysis of a generic braneworld model with induced 
gravity described by action (CQ). It was shown previously in [10|] that a low-density braneworld 
can mimic the expansion properties of the LCDM model. In particular, a universe consisting 
solely of baryons with fib — 0.04 can mimic the LCDM cosmology with a much larger 
'effective' value of the matter density fi^ CDM ~ 0.2—0.3. This property fixes the parameter 
(3 defined in © to be /3 « -5/4. 

The general analysis of the braneworld model with induced gravity performed in Sec. II III 
demonstrates that, in the neighborhood of the source, at distances smaller than r* given by 
( TT6l) . the solution is determined mainly by the quadratic part Q a b in (HI), which means that 
it respects the 'bare' Einstein equation m 2 G a b — T ab = to a high precision. This effect is 
sometimes described as the 'gravity filter' of the DGP model [3j, which screens the scalar 
graviton in the neighborhood of the source making the gravity effectively Einsteinian. 

At large distances, the modified gravitational field equations can be used to explain the 
properties of the galactic rotation curves. Using the smallness of the rotational velocity v c , 
we have developed a perturbative scheme for reconstructing the metric in this region. In the 
leading order of expansion, in the region given by (ED, our result reproduces that obtained 
for the Randall-Sundrum braneworld model in 6|, |7|. The vacuum gravity in this case can 
be approximately described by the Einstein equations with the traceless projection C a b of 
the bulk Weyl tensor to the brane playing the role of some effective stress-energy tensor of 
galactic fluid (jSJ). The left inequality in ( T4T1) is satisfied in a real spiral galaxy such as our 
Milky Way with the distance r ~ 3 kpc at which the rotation velocity becomes constant, 
u c ~7x 10~ 4 , if £ < + (3\ v~ l r ~ 2 Mpc. In this paper, we restricted ourselves to this 
case only. In the opposite case, one has to solve a nonlinear differential equation (1291) . 

The role of dark matter in the braneworld model is played by the traceless tensor of the 



Weyl fluid. Using the methods of 12] , in which it was proposed to use the existing data on 
rotation curves and lensing measurements to constrain the equation of state of galactic fluid, 
we have shown that deflection angles in our model are about 75% of the usual predictions 
based on dark matter. This result confirms the analysis done in jlj] and allows for a direct 
verification of the braneworld paradigm. 

In the region near to the gravitational object, we derive a closed system of equations for 
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a static spherically symmetric situation in the case of zero anisotropic stress in the tensor 
C a b. We demonstrate that the effective braneworld equations in this case have the form of 
Einstein equations with a special modification on the right-hand side [see (177j) -(J79l)]. We 
find the Schwarzschild metric to be an approximate vacuum solution of these equations at 
distances well below the scale r*, in agreement with t he g eneral consideration of Sec. IIHI 



Expression (fT6|) for generalizes the length scale 11] of the DGP model, below which 
nonlinear effects become important, to the case of nonzero brane tension (nonzero (3) and 
bulk cosmological constant satisfying the Randall-Sundrum constraint Ars = 0. For the 
parameter values (3 ~ —5/4 and t ^ 2 Mpc determined from the explanation of dark- matter 
phenomena, the value of r* for the Sun turns out to be of the order of 1 pc, which is 
quite large so that the solar-system experiments are in accord with the general-relativistic 
expectations. 

The density scale ( 1741) which appears in our model is estimated to be pl ~ 10 -22 g-cm~ 3 . 
This value lies far below the density of compact objects, such as stars and planets. Thus, 
inside such objects, we expect only slight modification of Newtonian dynamics due to the 
bulk effects. However, the local disk density of a typical spiral galaxy is pdisk ~ 3—12 x 
10 -24 g • cm -3 , and, therefore, the bulk effects inside the disk should be significant. The 
interior problem for such objects remains a challenge for future investigation. The solution 
to this problem will hopefully give a relation between the asymptotic rotational velocity v c 
and the mass of a spiral galaxy. It may help to verify the Tully-Fisher relation and extend 
our results to a wider class of astrophysical objects such as ellipticals and ultra compact 
dwarf galaxies. 
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